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Abstract 

We  compare  the  results  of  strain  state  analysis  (SSA)  and  photoluminescence  (PL) 
of  six  InGaN/GaN  quantum  well  samples  with  un-doped,  well-doped,  and 
barrier-doped  structures.  Based  on  the  SSA  images,  a  strain  relaxation  model  is 
proposed  for  describing  the  nanostructure  differences  between  the  three  sets  of 
sample  of  different  doping  conditions.  In  the  barrier-doped  samples,  the 
hetero-structure-induced  strains  are  fully  relaxed  such  that  spinodal  decomposition 
is  effectively  induced.  Therefore,  strongly  clustering  nanostructures  (quantum  dots) 
are  observed.  In  the  well-doped  samples,  strains  are  partially  relaxed  and  the 
spinodal  decomposition  process  can  be  slightly  induced.  Hence,  weaker 
composition  fluctuations  are  observed.  Then,  in  the  un-doped  samples,  the 
un-relaxed  strains  result  in  higher  miscibility  between  InN  and  GaN,  leading  to  the 
relatively  more  uniform  composition  distributions.  Between  the  low-  and 
high-indium  samples,  higher  indium  content  leads  to  a  stronger  clustering  behavior. 
The  strain  relaxations  in  the  well-doped  and  barrier-doped  samples  result  in  their 
unclear  S-shape  behaviors  of  PL  spectral  peaks.  The  enhanced  carrier  localization 
and  reduced  quantum-confined  Stark  effect  in  the  barrier-doped  samples  are 
responsible  for  their  significant  increases  of  radiative  efficiency. 
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I.  Introduction 

It  has  been  shown  that  different  doping  conditions,  including  different  doping 
layers,  in  InGaN/GaN  quantum-well  (QW)  structures  led  to  different  optical 
characteristics.  Such  effects  of  silicon  doping  have  been  widely  studied.  Several 
models  were  proposed  for  interpreting  the  optical  and  material  characteristics  of 
InGaN/GaN  QWs  with  silicon  doping.  The  most  commonly  used  model  is  the 
carrier  screening  effect  of  the  strain-induced  piezoelectric  field  and  hence  the 
reduction  of  the  quantum-confined  Stark  effect  (QCSE).  This  effect  leads  to  the 
reduced  Stokes  shift  [1-3],  the  decrease  of  radiative  recombination  lifetime 
[2-4,10],  the  blue  shift  of  PL  spectrum  [3,5],  and  the  enhancement  of  photon 
emission  efficiency  [1,3,6].  Also,  material  analyses  have  led  to  the  conclusions  of 
growth  mode  change  [6],  nanostructure  alteration  [7],  formation  of  abrupt  QW 
interfaces  [2,8],  strain  relaxation  [7,9],  and  higher  potential  uniformity  in  QWs  [3]. 
Meanwhile,  it  has  been  shown  that  different  barrier-doping  concentrations  in 
InGaN/GaN  QW  structures  resulted  in  different  optical  characteristics  [11].  These 
results  are  not  necessarily  mutually  consistent,  particularly  in  the  material 
nanostructures  of  such  silicon-doped  QW  samples. 

In  particular,  due  to  the  large  lattice  mismatch  between  InN  and  GaN, 
indium-rich  clusters  are  usually  formed  in  such  a  QW  structure  through  spinodal 
decomposition  and  other  effects  [12-14].  In  some  situations,  quantum-dot-like 
nanostructures  can  be  observed  [14-16].  With  different  silicon  doping  conditions, 
the  clustering  structures  are  expected  different.  Recently,  the  techniques  of 
high-resolution  transmission  electron  microscopy  (HRTEM)  have  been  widely 
used  for  understanding  the  properties  of  threading  dislocations  [17,18]  and 
nano-structures  in  GaN  and  InGaN/GaN  QW  samples.  In  particular,  the  existence 
of  indium-rich  clusters  in  InGaN/GaN  QW  structures  has  been  an  interesting  issue 
of  discussion  [19-22].  Although  some  researchers  claimed  that  the  observations  of 
clusters  were  due  to  the  artifacts  [19,20],  created  through  electron  beam  exposure 
during  HRTEM  measurements,  the  existence  of  clusters  really  relies  on  the 
conditions  of  sample  preparation.  Under  different  conditions  of  crystal  growth, 
clusters  may  or  may  not  exist.  Also,  the  cluster  characteristics  depend  on  those 
conditions. 

In  this  research,  we  use  the  strain-state  analysis  (SSA)  method  [23]  to  study 
the  nanostructures  of  InGaN/GaN  QWs  of  different  doping  conditions  (un-doped, 
well-doped,  and  barrier-doped)  and  different  average  indium  contents.  It  is  found 
that  the  clustering  phenomenon  is  the  strongest  in  the  barrier-doped  samples  and 
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the  weakest  in  the  un-doped  samples.  Among  the  samples  of  different  average 
indium  contents,  clustering  behavior  is  relatively  stronger  in  the  high-indium 
samples.  Based  on  the  SSA  results,  the  detailed  nanostructures  are  analyzed  and  the 
relative  indium  contents  can  be  confirmed.  Besides,  optical  properties  based  on 
photoluminescence  (PL)  measurements  are  provided  to  show  the  consistent  trends 
in  the  comparison  of  different  doping  conditions. 

II.  Sample  Descriptions  and  Experimental  Procedures 

The  six  InGaN/GaN  QW  samples  of  different  doping  conditions  and  different 
average  indium  contents  were  prepared  with  MOCVD  growth  on  sapphire  (0001) 
substrate.  They  all  consisted  of  five  QW  periods,  with  2.5  nm  in  well  width  and  7.5 
nm  in  barrier  width,  on  top  of  an  un-doped  GaN  buffer  layer  of  1.52  pm  in 
thickness.  The  growth  temperatures  were  1100  and  800  °C  for  the  GaN  barriers  and 
InGaN  wells,  respectively.  The  doping  concentration  of  silicon  was  5x10  cm' , 
either  in  barriers  or  wells.  They  are  classified  into  two  groups  of  low  and  high 
average  indium  contents.  The  low-indium  set  of  samples  has  an  estimated  content 
of  1 5  %.  The  high-indium  set  has  indium  content  of  25  %  in  estimation.  In  the 
low-indium  set,  the  samples  of  un-doped,  well-doped  and  barrier-doped  conditions 
were  denoted  with  samples  LU,  LW,  and  LB,  respectively.  The  counterparts  of  the 
high-indium  set  are  assigned  as  samples  HU,  HW,  and  HB,  respectively. 

The  HRTEM  investigations  (for  SSA)  were  performed  using  a  Philips  Tecnai 
F30  field-emission  electron  microscope  using  an  accelerating  voltage  of  300  kV 
and  a  probe  forming  lens  of  Cs  =  1.2  mm.  All  the  high-resolution  micrographs 
were  taken  with  the  two-electron-beam  interference.  The  current  density  focused 
onto  the  samples  was  estimated  to  be  lower  than  16  A/cm  .  With  the  two-beam 
interference  data,  we  could  conduct  the  SSA  for  composition  distribution  images 
[23].  In  the  HRTEM  operation,  to  avoid  the  possible  influence  of  electron 
illumination  on  the  recorded  results  [19,20],  the  two-beam  interference  images 
were  always  taken  within  5-10  seconds  after  the  electron  beam  was  focused  onto 
the  target  spots.  Our  experiments  showed  that  the  processed  SSA  images  were  not 
significantly  changed  within  30  seconds  of  the  illumination  of  focused  electron 
beams.  To  obtain  the  SSA  images,  the  [0001]  lattice  size  distributions  were  first 
calibrated  from  the  lattice  fringe  images.  Then,  based  on  Vegard’s  law  for  the 
relation  between  the  lattice  size  and  indium  composition,  composition  images 
could  be  obtained.  These  compositions  were  estimated  based  on  the  assumption  of 
a  specimen  thickness  larger  than  50  nm.  PL  measurements  were  carried  out  with 
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the  325  nm  line  of  a  35  mW  He-Cd  laser.  The  samples  were  placed  in  a  cryostat  for 
temperature-dependent  measurements  ranging  from  10  to  300  K. 

III.  Images  of  Strain  State  Analysis 

Figs.  1-6  show  the  typical  SSA  results  of  samples  LU,  LW,  LB,  HU,  HW,  and 
HB,  respectively.  In  each  figure,  part  (a)  shows  the  lattice  fringe  image  of 
two-beam  interference  along  the  [0001]  direction  and  part  (b)  shows  the  processed 
SSA  image.  One  can  see  that  the  lattice  fringe  image  cannot  clearly  describe  the 
nanostructures  of  materials.  The  colors  in  the  SSA  image  represent  the  estimated 
indium  composition  distributions.  The  relationship  between  the  lattice  parameter 
ratio  dj/doaN  (dGaN  and  dj  for  the  lattice  sizes  of  GaN  and  InGaN,  respectively)  and 
the  indium  composition  Xjn  (%)  is  shown  in  the  legend  at  the  bottom  of  the  figure. 
It  is  noted  that  their  relationship  is  not  linear,  particularly  when  the  indium 
composition  is  high.  In  each  SSA  image,  line  scan  profiles  of  indium  composition 
along  the  indicated  white  lines  are  plotted.  The  scales  for  the  profiles  in  the  two 
dimensions  are  the  same.  In  each  SSA  image,  a  QW  surrounded  by  the  two  barriers 
are  illustrated. 

Fig.  1(b)  shows  the  SSA  image  of  sample  LU.  One  can  see  that  the  QW  layer 
(green  region)  is  quite  uniform.  The  indium  clusters  (reddish  regions)  within  the 
QW  layer  are  quite  small  and  shallow.  The  GaN  barriers  (deep  blue  regions)  are 
clear.  It  means  that  indium  is  well  confined  within  the  well.  From  the  horizontal 
line-scan  profile,  one  can  observe  quite  a  weak  fluctuation  in  indium  composition 
along  the  well  layer.  The  fluctuation  contrast  (the  difference  between  the  maximum 
and  the  minimum  in  the  scan  range)  is  in  the  range  of  0.03-0.04  (18-24  %)  along 
the  QW  layer.  In  Fig.  2(b)  for  sample  LW,  we  can  observe  that  more  indium 
clusters  are  formed  in  the  QW  layer.  Also,  more  indium  atoms  diffuse  from  the 
well  into  barriers.  Here,  the  fluctuations  of  the  line-scan  profiles  are  relatively 
stronger  in  comparing  with  those  of  sample  LU.  The  fluctuation  contrast  is  now  in 
the  range  of  0.04-0.05  (22-28  %).  In  Fig.  3(b)  for  sample  LB,  we  can  see  two  clear 
indium-rich  clusters  in  the  QW  (two  large  red  spots).  From  the  line-scan  profiles, 
one  can  see  that  the  indium  composition  fluctuation  in  sample  LB  is  much  stronger 
than  those  of  the  other  two  samples.  The  fluctuation  contrast  is  now  around  0.08 
(44  %). 

As  shown  in  Fig.  4(b)  for  sample  HU,  although  the  QW  interface  is  blurred, 
indium  is  basically  confined  within  the  well.  Here,  within  the  well  a  red  spot  of 
indium  aggregation  can  be  observed.  Similar  to  sample  LU,  one  can  also  observe  a 
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relatively  weaker  fluctuation  in  indium  composition  along  the  well  layer  in  the 
line-scan  profile.  The  fluctuation  contrast  is  about  0.04  (22  %)  along  the  QW  layer. 
In  Fig.  5(b)  for  the  SSA  image  of  sample  HW,  the  QW  is  not  as  well  shaped  as 
sample  HU.  The  indium  composition  fluctuations  in  both  directions  are  relatively 
stronger  in  comparing  with  sample  HU.  In  particular,  more  indium-aggregated 
clusters  can  be  observed  within  the  well  layer.  The  fluctuation  contrast  is  also 
around  0.04  (22  %)  along  the  QW  layer.  As  shown  in  Fig.  6(b),  the  SSA  image  of 
sample  HB  shows  quite  a  different  nanostructure  from  the  other  five  samples.  Here, 
the  QW  layer  becomes  unclear.  Instead,  a  distribution  of  clusters  of  different  sizes 
and  shapes  exists.  The  indium  composition  fluctuation  is  much  stronger  than  those 
of  samples  HU  and  HW.  The  contrast  is  now  around  0.08  (44  %)  along  the  QW 
layer. 

A  question  has  been  raised  about  the  possible  influence  of  electron  beam 
exposure  on  the  material  nanostructures  during  the  measurement  of  transmission 
electron  microscopy  [19,20].  We  have  tested  such  a  possibility  by  varying  the 
electron  exposure  time.  The  SSA  images  of  sample  LB  shown  in  Figs.  7(a)  and  (b) 
were  obtained  under  electron  beam  exposure  for  less  than  10  sec  and  about  two 
minutes,  respectively.  Fig.  7(a)  is  different  from  Fig.  3(b)  because  they  were 
obtained  at  different  locations  of  sample  LB.  Between  Figs.  7(a)  and  (b),  one  can 
see  that  although  the  details  are  slightly  changed,  the  basic  structures  are  the  same. 
In  other  words,  with  up  to  two-minute  exposure,  we  can  still  obtain  the  required 
information  about  the  nanostructures  of  the  samples.  Normally,  when  we  operated 
the  HRTEM  system,  it  took  less  than  10  sec  after  the  electron  beam  was  focused 
onto  a  target  spot  to  obtain  an  image.  In  our  research,  we  have  done  our  best  to 
reduce  the  electron  beam  effect.  Hence,  the  SSA  images  in  this  research  can  truly 
represent  the  sample  nanostructures. 

In  Figs.  8  and  9,  we  show  the  average  line-scan  profiles  in  the  [0001]  direction, 
averaged  over  the  image  range  along  the  [11-20]  direction,  plotted  as  functions  of 
distance  for  the  low-  and  high-indium  samples,  respectively.  Two  to  three  SSA 
images  like  Figs.  1  -6  at  different  locations  in  each  sample  were  used  for  taking  the 
average.  One  can  see  that  generally  the  well  widths  of  the  high-indium  samples  are 
larger  than  those  of  the  low-indium  ones.  One  can  also  observe  that  the  interfaces 
between  the  well  and  barriers  of  the  low-indium  samples  are  generally  shaper  than 
those  of  the  high-indium  ones.  Meanwhile,  generally  the  indium  concentrations  at 
the  peaks  of  the  high-indium  samples  are  higher  than  those  of  the  low-indium  ones. 
In  addition,  in  the  low-indium  samples,  more  indium  atoms  diffuse  from  the  well 
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into  the  lower  barrier.  Indium  out-diffusion  in  the  high-indium  samples  is  more 
symmetric. 

The  average  indium  content  of  each  sample  can  be  calibrated  by  integrating 
the  average  line-scan  results  in  Fig.  8  or  9  over  the  whole  range  and  then  dividing 
the  integrated  result  by  the  designated  QW  width  (2.5  nm).  The  whole-range 
integration  implies  the  coverage  of  the  indium  atoms,  which  diffuse  from  the  well 
into  barriers.  Such  results  of  average  indium  content  become  more  accurate  if  more 
SSA  images  like  Figs.  1-6  are  used  for  average.  Although  the  current  results  cannot 
be  very  accurate,  they  can  still  be  used  for  a  general  comparison.  The  calibrated 
average  indium  contents  of  samples  LU,  LW,  and  LB  are  14.3,  17.1,  and  19.2  %, 
respectively.  Those  of  samples  HU,  HW,  and  HB  are  25.6,  26.7,  and  26.9  %, 
respectively.  These  results  are  generally  consistent  with  the  designated  indium 
contents  of  15  and  25  %  for  the  low-  and  high-indium  samples,  respectively.  The 
generally  higher  indium  contents  from  calibration  for  the  barrier-  and  well-doped 
samples  can  be  attributed  to  the  stronger  indium  diffusion  into  barriers  in  these  two 
kinds  of  sample.  With  more  indium  atoms  in  the  barriers,  the  calibration  for  indium 
composition  in  processing  the  SSA  images  can  be  overestimated.  However,  the 
possibility  of  really  higher  average  indium  content  in  the  doped  samples  cannot  be 
ruled  out.  With  different  strain  conditions  (due  to  different  doping  conditions),  the 
incorporation  efficiency  of  indium  atom  can  be  different  even  the  flow  rate  is  the 
same  during  growth. 

IV.  Optical  Characterization  Results 

To  understand  the  corresponding  optical  properties  of  these  samples,  basic 
optical  characterizations  were  performed.  Fig.  10  shows  the  PL  spectral  peak 
positions  as  functions  of  temperature  for  all  the  six  samples  with  the  filled  and 
empty  symbols  for  the  low-  and  high-indium  samples,  respectively.  Clear  S-shape 
variations  can  be  seen  in  the  curves  for  samples  LU  and  HU.  On  the  contrary,  the 
S-shape  behaviors  are  unclear  in  the  well-  and  barrier-doped  samples.  The  PL  peak 
positions  of  the  un-doped  and  well-doped  samples  in  either  group  of  sample  are 
quite  close,  particularly  in  the  low-temperature  range.  The  significant  blue  shifts  of 
the  PL  peaks  of  samples  LB  and  HB,  with  respect  to  those  of  the  un-doped  and 
well-doped  samples  show  a  clear  feature  of  barrier  doping. 

Fig.  11  shows  the  integrated  PL  intensities  as  functions  of  temperature  of 
samples  LU,  LW,  and  LB.  Fig.  12  shows  the  counterparts  of  samples  HU,  HW,  and 
HW.  The  thermal  quenching  behavior  in  integrated  luminescence  intensity 
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represents  the  radiative  efficiency.  From  Figs.  11  and  12,  one  can  see  that  the 
barrier-doped  samples  always  have  the  highest  radiative  efficiency,  followed  by  the 
well-doped  samples  and  then  the  un-doped  samples.  Hence,  in  our  samples,  silicon 
doping,  particularly  doping  in  barriers,  can  lead  to  a  higher  efficiency  of 
spontaneous  emission. 

V.  Discussions 

The  SSA  images  show  that  barrier  doping  results  in  stronger  indium  clustering. 
This  result  may  be  caused  by  the  significant  relaxation  of  hetero-structure-induced 
strain  (the  strain  in  QW  due  to  the  lattice  mismatch  between  wells  and  barriers)  in 
the  barrier-doped  samples,  particularly  in  sample  HB,  through  certain  mechanisms. 
In  sample  HB,  with  the  growth  temperature  at  800  °C  and  nominal  indium  content 
about  25  %,  the  sample  condition  falls  into  the  regime  of  spinodal  decomposition 
(the  unstable  regime)  in  the  T-x  phase  diagram  of  ternary  InxGai.xN  compounds  for 
relaxed  and  strained  layers  with  the  interface  orientation  perpendicular  to  the 
hexagonal  axis  of  the  crystal  [12,24].  In  this  regime,  spinodal  decomposition  and 
phase  separation  can  be  initiated  by  a  small  local  compositional  fluctuation.  Hence, 
strongly  clustering  behaviors  can  be  observed  in  sample  HB.  The  mechanisms  for 
relaxing  the  hetero-structure-induced  strain  can  be  the  generation  of  local 
dislocations,  which  is  caused  by  the  enhanced  lattice  mismatch  between  a  well  and 
the  surrounding  barriers  because  of  the  substitutions  of  a  high-concentration  of 
larger  gallium  atoms  (atomic  weight  as  31)  by  smaller  silicon  atoms  (atomic 
weight  as  14)  in  the  barriers.  The  HRTEM  images  in  the  similar  samples  of  our  and 
other’s  previous  research  have  shown  the  existence  of  local  dislocations  in  an 
InGaN/GaN  QW  sample  of  a  high  nominal  indium  content  [25,26].  The  other 
mechanism  for  strain  relaxation  in  the  barrier-doped  sample  is  the  existence  of 
plenty  carriers  in  the  barriers  and  wells  (carriers  flow  into  wells)  due  to  the  heavy 
silicon  doping  in  the  barriers.  The  carriers  screen  the  piezoelectric  fields  and  affect 
the  strain  energy  distributions.  However,  it  is  believed  that  the  major  mechanism 
for  strain  relaxation  is  the  formation  of  local  dislocations.  The  discussions 
described  above  fit  well  the  conditions  of  sample  HB.  In  sample  LB,  because  the 
nominal  indium  content  is  lower  (about  15  %),  the  hetero-structure-induced  strain 
is  not  as  strong  as  that  in  sample  HB.  Therefore,  the  condition  may  fall  into  a 
region  between  the  unstable  and  metastable  regimes  in  the  T-x  phase  diagram.  In 
this  situation,  spinodal  decomposition  and  phase  separation  processes  can  still  be 
induced.  However,  these  processes  in  the  meta-stable  regime  take  a  longer  time 
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than  that  in  the  unstable  regime,  leading  to  a  relatively  weaker  clustering  behavior. 

In  the  well-doped  samples,  because  of  the  substitutions  of  larger  gallium  or 
indium  atoms  by  smaller  silicon  atoms  (lattice  mismatch  reduced)  and  the  effect  of 
carrier  screening,  the  hetero-structure-induced  strain  has  been  partly  relaxed  such 
that  the  conditions  may  belong  to  the  meta-stable  regime  in  the  T-x  phase  diagram. 
Hence,  certain  clusters  were  generated  although  the  clustering  behaviors  are  not  as 
strong  as  those  in  the  barrier-doped  samples.  In  the  un-doped  samples,  the  strong 
hetero-structure-induced  strains  are  preserved  such  that  the  conditions  are 
essentially  located  in  the  stable  regime.  In  this  situation,  the  solubility  of  InN  in 
GaN  is  higher  such  that  the  clustering  behaviors  are  the  weakest  among  the  three 
kinds  of  sample.  It  is  noted  that  with  a  higher  nominal  indium  content,  the 
condition  in  the  phase  diagram  can  be  closer  to  the  meta-stable  or  even  spinodal 
decomposition  regime  and  hence  the  clustering  behavior  can  be  stronger. 

From  the  SSA  results  above,  we  can  build  a  model  for  potential  variation 
along  a  quantum  well  layer.  In  the  situations  of  samples  HU  and  LU,  potential 
fluctuations  are  weak.  A  potential  minimum  is  surrounded  by  secondary  minima  of 
shallow  barriers  in  between.  In  the  cases  of  samples  HW  and  LW,  potential 
variations  are  relatively  stronger.  In  this  situation,  a  potential  minimum  is  also 
surrounded  by  secondary  minima.  However,  the  barriers  between  the  minima  are 
now  relatively  higher  such  that  certain  energy  is  required  for  carrier  transport 
between  the  minima.  Then,  for  samples  HB  and  LB,  sharp  minima  with  few 
secondary  dips  are  assumed.  Therefore,  carriers  can  relax  down  to  the  absolute 
potential  minima  without  a  cascading  process  through  the  secondary  minima.  The 
difference  between  the  well-doped  and  un-doped  samples  is  that  more  energy  is 
required  for  overcoming  the  barriers  in  samples  LW  and  HW. 

The  relative  larger  PL  photon  energies  in  the  barrier-doped  samples  can  be  due 
to  the  reduction  of  piezoelectric  fields  that  originates  from  strain  relaxation  and 
carrier  screening.  Also,  it  can  be  due  to  the  increase  of  quantum  confinement, 
which  results  from  the  formation  of  nano-scale  clusters.  In  the  well-doped  samples, 
more  carriers  in  the  well  layers  should  result  in  significant  increases  of  PL  photon 
energies  due  to  strong  carrier  screening.  However,  the  PL  photon  energies  of  the 
well-doped  samples  are  quite  close  to  those  of  the  un-doped  samples.  Such  a  result 
can  be  attributed  to  the  ineffective  carrier  screening  effect.  It  can  also  be 
interpreted  as  the  result  of  the  counteraction  between  band  gap  renormalization 
(red-shift  trend)  and  piezoelectric  field  reduction  (blue-shift  trend).  As  shown  in 
Fig.  10,  the  temperature-dependent  PL  spectral  peak  energies  reveal  clear  S-shape 
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behaviors  in  samples  LU  and  HU.  Such  an  S-shape  behavior  is  usually  attributed  to 
carrier  localization  [27,28]  and/or  QCSE  [29,30].  Because  the  S-shape  behavior  is 
unclear  in  either  well-doped  or  barrier-doped  samples,  whose  strains  are  supposed 
to  be  partially  or  fully  relaxed,  the  major  mechanism  for  the  S-shape  behavior  in 
our  samples  must  be  the  strain  distribution  or  QCSE.  Without  significant  relaxation 
of  strains  in  the  un-doped  samples,  the  S-shape  behaviors  are  quite  prominent.  The 
significant  enhancement  of  radiative  recombination  efficiency  in  the  barrier-doped 
samples  can  be  attributed  to  the  strong  carrier  localization  (localized  from 
non-radiative  recombination  centers)  and  the  reduced  QCSE  (improved  wave 
function  overlap  between  electron  and  hole). 

VI.  Conclusions 

In  summary,  we  have  compared  the  results  of  SSA  and  PL  of  six  InGaN/GaN 
QW  samples  with  un-doped,  well-doped,  and  barrier-doped  structures.  Based  on 
the  SSA  images,  a  strain  relaxation  model  was  built  for  describing  the 
nanostructure  differences  between  the  three  sets  of  sample  of  different  doping 
conditions.  In  the  barrier-doped  samples,  the  hetero-structure-induced  strains  were 
fully  relaxed  such  that  spinodal  decomposition  was  effectively  induced.  Therefore, 
strongly  clustering  nanostructures  were  observed.  In  the  well-doped  samples, 
strains  were  partially  relaxed  and  the  spinodal  decomposition  process  could  be 
slightly  induced.  Hence,  weaker  composition  fluctuations  were  observed.  Then,  in 
the  un-doped  samples,  the  un-relaxed  strains  resulted  in  higher  miscibility  between 
InN  and  GaN,  leading  to  the  relatively  more  uniform  composition  distributions. 
Between  the  low-  and  high-indium  samples,  higher  indium  content  led  to  a 
stronger  clustering  behavior.  The  strain  relaxations  in  the  well-doped  and 
barrier-doped  samples  resulted  in  their  unclear  S-shape  behaviors  of  PL  spectral 
peaks.  The  enhanced  carrier  localization  and  reduced  QCSE  in  the  barrier-doped 
samples  were  responsible  for  their  significant  increases  of  radiative  efficiency. 
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Fig.  1  Typical  HRTEM  lattice  fringe  image  (a)  and  the  SSA  image  (b)  after 

software  process  of  samples  LU. 
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Fig.  2  Typical  HRTEM  lattice  fringe  image  (a)  and  the  SSA  image  (b)  after 

software  process  of  samples  LW. 
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Fig.  3  Typical  HRTEM  lattice  fringe  image  (a)  and  the  SSA  image  (b)  after 

software  process  of  samples  LB. 
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Fig.  4  Typical  HRTEM  lattice  fringe  image  (a)  and  the  SSA  image  (b)  after 

software  process  of  samples  HU. 
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Fig.  5  Typical  HRTEM  lattice  fringe  image  (a)  and  the  SSA  image  (b)  after 

software  process  of  samples  HW. 
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Fig.  7  SSA  images  of  sample  LB  under  different  electron-beam  exposure  times:  (a) 

<10  sec  and  (b)  120  sec. 
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Fig.  8  Average  indium  composition  profiles  along  the  crystal  growth  direction  for 

samples  LU,  LW,  and  LB. 
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Fig.  9  Average  indium  composition  profiles  along  the  crystal  growth  direction  for 

samples  HU,  HW,  and  HB. 


22 


Temperature  (K) 


Fig.  10  PL  spectral  peak  variations  with  temperature  of  the  six  samples. 
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Fig.  11  Integrated  PL  intensity  variations  with  temperature  of  samples  LU,  LW,  and 

LB. 
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Fig.  12  Integrated  PL  intensity  variations  with  temperature  of  samples  HU,  HW, 

and  HB. 
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Temperature-dependent  femtosecond  pump-probe  experiments  are  performed  to  explore  the 
ultrafast  carrier-relaxation  processes  in  an  InGaN/GaN  quantum-well  sample,  in  which 
nanometer-scale  cluster  structures  have  been  identified.  Combined  with  the  time-resolved 
photoluminescence  results,  we  can  identify  three  stages  of  carrier  relaxation.  The  fast-decay  time, 
ranging  from  several  hundred  femtoseconds  to  1  picosecond,  corresponds  to  the  process  reaching  a 
local  quasi-equilibrium  condition,  in  which  carriers  reach  a  thermal  distribution  within  one  or  a  few 
nearby  indium-rich  clusters.  The  slow-decay  time,  ranging  from  tens  to  a  couple  hundred 
picoseconds,  corresponds  to  the  process  reaching  a  global  quasi-equilibrium  condition,  in  which 
carriers  reach  a  thermal  distribution  among  different  clusters  of  various  potential  minima.  In  this 
stage,  the  mechanism  of  carrier  transport  over  barriers  between  clusters  dominates  the  relaxation 
process.  Finally,  carrier  recombination  dominates  the  relaxation  process  with  the  carrier  lifetime  in 
the  range  of  a  few  nanoseconds.  ©  2004  American  Institute  of  Physics,  [DOI:  10.1063/1.1784033] 


InGaN/GaN  quantum- well  (QW)  structures  have  been 
widely  used  for  blue-green  light-emitting  devices.  Due  to  the 
large  lattice  mismatch  between  InN  and  GaN,  indium-rich 
aggregation  usually  occurs  in  InGaN,  leading  to  the  forma¬ 
tion  of  nanometer-(nm)-scale  clusters  and  hence  potential 
fluctuations.1”6  Carriers,  once  generated,  are  trapped  in  the 
localized  states  of  potential  minima  and  recombine  effec¬ 
tively.  The  cluster  structure  is  the  key  to  the  efficient  emis¬ 
sion  in  such  a  compound.  Normally,  cluster  structures  with 
well-defined  potential  minima,  such  as  those  in  a  sample  of 
barrier  doping  with  silicon,  result  in  higher  emission 
efficiency.7  Because  the  information  of  carrier  relaxation  in 
such  a  sample  can  help  us  in  understanding  better  its  nano¬ 
structure  and  hence  the  recombination  process,  an  ultrafast 
pump-probe  study  is  useful  for  the  development  of  related 
devices.  Although  a  temperature-dependent  femtosecond  (fs) 
pump-probe  experiment  on  InGaN/ InGaN  QWs  of  low  in¬ 
dium  contents  was  reported,8  no  cluster-related  phenomena 
were  mentioned. 

In  this  letter,  we  report  the  results  of  temperature- 
dependent  fs  pump-probe  experiments  on  an  InGaN/GaN 
QW  sample.  By  combining  the  results  of  picosecond-scale 
time-resolved  photoluminescence  (TRPL)  measurement,  we 
can  identify  different  time  scales  of  carrier  relaxation  pro¬ 
cess.  The  sample  consists  of  five  periods  of  QW  with  2.5  nm 
In0 15Ga0  g5N  wells  and  10  nm  GaN  barriers.  The  details  of 
sample  structure  and  growth  condition  can  be  found  in  Ref. 
9.  Cluster  structures  in  the  nm  scale  have  been  identified  in 
this  sample  with  high-resolution  transmission  electron 
microscopy.10  In  the  pump-probe  experiment,  a  76  MHz, 
100  fsTirsapphire  laser,  pumped  with  a  Nd:  YV04  laser,  was 
frequency  doubled  through  a  0.5-mm-thick  BBO  crystal.  The 
second-harmonic  pulses  were  split  into  two  beams  for  degen¬ 
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erate,  noncollinear,  and  cross-polarized  pump-probe  mea¬ 
surements.  A  half-wave  plate  and  two  polarizers  were  used  to 
minimize  coherent  artifacts.  The  1  jim  stepping  resolution  of 
the  computer-controlled  translation  delay  stage  accounts  for 
7  fs  step  size  in  the  time  domain.  Pump  photon  energies  at 
2.98  and  3.14  eV  were  chosen,  as  explained  in  the  following* 
with  average  powers  at  12  and  32  mW  (4.96  and 
13.2  ^tJ/cm2),  respectively.  The  sample  was  back-side  pol¬ 
ished  to  guarantee  effective  transmission  of  the  probe  beam. 
Figure  1  shows  the  photoluminescence  (PL)  and  photolumi¬ 
nescence  excitation  (PLE)  spectra  at  10  K.  At  around 
3.14  eV,  there  is  an  absorption  peak  in  the  PLE  spectrum. 
Also,  near  2.98  eV,  the  PLE  intensity  reaches  its  minimum 
before  it  rapidly  increases  as  the  photon  energy  falls  into  the 
PL  spectral  range.  We  choose  these  two  photon  energies  for 
pumping  the  sample  (see  the  arrows  in  Fig.  1). 

In  Fig.  2,  we  show  the  differential  transmission,  1ST  IT, 
as  functions  of  time  at  various  sample  temperatures  when  the 


FIG.  1.  PL  and  PLE  spectra  at  10  K.  The  used  pump  photon  energies  are 
indicated  with  arrows. 
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FIG.  2,  Time  evolution  of  differential  transmission  at  various  sample  tem¬ 
peratures  when  the  pump  energy  is  2.98  eV. 

pump-probe  photon  energy  is  2.98  eV.  One  can  see  that  in 
each  temporal  trace,  after  the  pump  pulse,  differential  trans¬ 
mission  decays  with  time  first  quite  fast  in  the  fs  range,  fol¬ 
lowed  by  a  slower  picosecond-  (ps)-scale  decay.  Traces  in 
Fig.  2  show  the  temperature-dependent  variations  of  the  fast- 
decay  speed  and  relative  level  of  the  slow-decay  range.  The 
decay  behavior  around  105  K  is  quite  different  from  those  at 
lower  or  higher  temperatures.  In  particular,  oscillations  in  the 
slow-decay  range  can  be  clearly  seen.  Such  oscillations  are 
supposed  to  be  due  to  the  existence  of  coherent  acoustic 
phonons.11,12 

To  explore  the  decay  phenomena,  we  evaluated  the  ratio 
of  differential  transmission  of  the  level  in  the  slow-decay 
range  (we  used  the  value  at  4  ps  to  represent  this  level)  over 
that  of  the  peak.  The  results  are  shown  in  the  curve  with 
closed  triangles  in  Fig.  3.  The  vertical  line  segments  repre¬ 
sent  the  error  bars,  which  were  obtained  from  ten  times  mea¬ 
surements.  Also  shown  in  Fig.  3  is  the  curve  with  closed 
circles  for  the  PL  spectral  peak  position  as  a  function  of 
temperature.  The  PL  spectral  peak  position  shows  a  typical 
S-shape  variation,  which  is  a  key  feature  of  a  clustering 
structure,  indicating  the  effect  of  carrier  localization.9,13,14  In 


Temperature  (K) 

FIG.  3.  Temperature-dependent  variations  of  PL  spectral  peak  and  transmis¬ 
sion  ratios  at  two  pump  energies. 

Downloaded  18  Aug  2004  to  128.220.144.249.  Redistribution  subject 


Wang  et  al. 


Temperature  (K) 

FIG.  4.  Temperature-dependent  variations  of  the  fast-  and  slow-decay  times 
and  the  carrier  lifetime  when  the  pump  energy  is  2.98  eV. 

Fig.  3,  it  is  interesting  to  observe  that  when  the  sample  is 
pumped  with  photon  energy  of  2.98  eV,  a  peak  of  transmis¬ 
sion  Tatio  exists  around  105  K,  at  which  the  PL  spectral  peak 
starts  to  show  the  significant  blueshift  trend.  Based  on  the 
observation,  we  can  build  a  model  as  follows:  At  low  tem¬ 
peratures  (<105  K),  the  pump  photon  energy  is  far  above 
the  peak  level  of  carrier  thermal  distribution.  In  this  situa¬ 
tion,  once  carriers  are  generated,  in  the  process  of  approach¬ 
ing  the  quasi-equilibrium  condition,  carriers  move  rapidly 
into  the  states  of  lower  levels.  Therefore,  after  the  fast  decay, 
the  carrier  density  at  the  pump  energy  level  is  low.  When 
temperature  increases  to  a  certain  level  (around  105  K  in  this 
sample),  the  peak  level  of  carrier  thermal  distribution  is 
raised  to  become  close  to  the  pump  energy  level.  In  this 
situation,  plenty  of  carriers  stay  in  the  levels  close  to  the 
pump  energy  during  the  relaxation  process  of  approaching  a 
quasi-equilibrium  condition.  Therefore,  the  transmission  ra¬ 
tio  shown  in  Fig.  3  reaches  a  maximum  around  105  K.  Be¬ 
yond  this  temperature,  the  phonon  effect  (leading  to  band 
gap  shrinkage  in  increasing  temperature)  gradually  domi¬ 
nates  and  the  peak  level  of  carrier  distribution  moves  below 
the  pump  energy  level.  Hence,  the  transmission  ratio  essen¬ 
tially  decreases  with  temperature. 

To  obtain  the  carrier  decay  times,  we  used  a  model  of 
double-exponential  variation,  convolved  with  a  Gaussian 
function  (the  pump  pulse),  to  fit  the  measured  data,15  as  il¬ 
lustrated  by  the  dashed  curves  in  Fig.  2.  With  such  fitting,  the 
fast-  and  slow-decay  times  were  obtained,  as  plotted  in  Fig. 
4.  One  can  see  that  the  fast  decay  times  (closed  circles)  are  in 
the  range  from  700  fs  to  1  ps.  It  generally  decreases  with 
temperature  with  a  shallow  maximum  around  105  K.  On  the 
other  hand,  the  slow-decay  times  (open  circles)  range  from 
50  to  200  ps  with  a  general  increasing  trend  with  tempera¬ 
ture.  A  minimum  around  105  K  can  be  clearly  seen.  In  Fig. 
4,  we  also  show  the  carrier  lifetimes,  which  were  obtained 
with  TRPL  measurement  by  using  a  streak  camera  of  about 
4  ps  in  temporal  resolution.  They  were  calibrated  from  the 
nanosecond-  (ns)-scale  decay  of  PL  intensity  when  the 
sample  was  excited  with  100  fs  pulses  at  3.14  eV.  These 
decay  times  are  in  the  range  of  a  few  ns  and  correspond  to 
the  carrier  lifetime  in  the  recombination  process.  In  the 
TRPL  measurement,  a  faster  decay  with  decay  time  in  the  ps 
range  was  also  observed.  This  ps-scale  decay,  corresponding 
to  the  open-circle  data  in  Fig.  4,  has  been  attributed  to  the 
carrier  transport  process  among  potential  minima  of  different 

levels.16  With  such  observations,  we  can  speculate  that  the 
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FIG.  5.  Temperature-dependent  variations  of  the  fast-  and  slow-decay  times 
when  the  pump  energy  is  3.14  cV. 

fast-decay  times  in  Fig.  4  describe  the  carrier  relaxation  pro¬ 
cess  toward  a  quasi-equilibrium  condition  in  a  localized  re¬ 
gion  (within  one  or  a  few  nearby  clusters).  Because  the 
sample  has  a  structure  between  two  and  zero  dimensions,  a 
quasi-continuous  state  distribution  still  exists.  Hence,  carrier 
redistribution  among  different  states  for  a  local  quasi¬ 
equilibrium  condition  dominates  carrier  dynamics  in  this 
time  range.  Such  a  process  becomes  faster  as  temperature 
increases,  as  clearly  described  in  Fig.  4.  After  this  stage, 
carriers  significantly  transport  to  neighboring  clusters  of 
higher  or  lower  potential  minima,  reaching  a  global  quasi- 
equilibrium  condition.  This  process  takes  an  average  time  of 
about  100  ps,  which  is  consistent  with  our  previous  calibra¬ 
tions  based  on  the  TRPL  measurement.9  At  the  temperature 
(around  105  K)  where  the  PL  spectral  peak  shows  a  signifi¬ 
cant  blue  shift,  such  a  relaxation  process  becomes  faster  be¬ 
cause  carriers  are  supposed  to  move  into  the  energy  states 
close  to  the  pump  level.  In  this  transport  process,  carriers 
need  to  gain  certain  thermal  energy  for  overcoming  the  bar¬ 
riers  between  neighboring  clusters.  After  the  transport  stage, 
electrons  with  the  global  quasi-equilibrium  condition  will 
eventually  relax  back  to  the  valence  band  with  a  time  scale 
of  carrier  lifetime. 

With  pump  at  3.14  eV,  the  transmission  ratios  are  shown 
as  the  curve  with  closed  squares  in  Fig.  3.  In  this  situation, 
the  peak  around  105  K,  like  that  in  the  case  of  2.98  eV 
pump,  disappears.  The  result  is  attributed  to  the  fact  that  the 
pump  energy  is  far  above  the  peak  level  of  thermal  carrier 
distribution,  even  with  the  blueshift  trend  of  PL  spectral  peak 
beyond  105  K.  The  fitted  fast-  and  slow-decay  times  in  the 
case  of  3.14  eV  pump  are  shown  in  Fig.  5.  The  results  are 
similar  to  those  in  Fig.  4  although  the  decay  times  with 
3,14  eV  pump  are  slightly  shorter.  It  has  been  pointed  out 
that  the  carrier  transport  process  through  overcoming  a  vari¬ 


ety  of  barrier  heights  should  follow  a  power-law  behavior  in 
time  evolution.17  However,  because  our  measurement  time 
range  (—10  ps)  represents  only  a  small  portion  of  the  whole 
slow-decay  process,  it  is  difficult  to  exactly  calibrate  its  time 
dependence. 

In  summary,  we  have  performed  temperature-dependent, 
fs  pump-probe  experiments  to  explore  the  ultrafast  carrier- 
relaxation  processes  in  an  InGaN/GaN  QW  sample,  in  which 
nm-scale  cluster  structures  have  been  identified.  Combined 
with  the  TRPL  results,  we  could  identify  three  stages  of  car¬ 
rier  relaxation.  The  fast-decay  time,  ranging  from  several 
hundred  fs  to  1  ps,  corresponded  to  the  process  toward  a 
local  quasi-equilibrium  condition,  in  which  carriers  reached 
a  thermal  distribution  within  one  or  a  few  nearby  clusters. 
The  slow-decay  time,  ranging  from  tens  to  a  couple  of  hun¬ 
dred  ps,  corresponded  to  the  process  toward  a  global  quasi¬ 
equilibrium  condition  Finally,  the  mechanism  of  carrier  re¬ 
combination  dominated  the  process  with  the  carrier  lifetime 
in  the  range  of  a  few  ns. 
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We  perform  femtosecond  degenerate  pump-probe  experiments  on  an  InGaN  thin  film  of  800  nm  in 
thickness.  The  observed  temperature-,  pump-photon-energy-,  and  pump-intensity-dependent 
variations  of  ultrafast  carrier  dynamics  manifest  the  variation  of  the  space-averaged  density  of  state 
with  energy  level  in  this  sample.  The  carrier  dynamics  is  controlled  by  the  shift  of  effective  band  gap 
and  hence  the  behavior  of  band  filling,  which  are  determined  by  the  combined  effect  of  band-gap 
renormalization  and  phonon  effect  (band-gap  shrinkage  with  increasing  temperature).  Two-photon 
absorption  and  free-carrier  absorption  can  be  observed  when  the  corresponding  density  of  state  is 
low  and  hence  the  band-filling  effect  is  weak.  The  variation  of  the  space-averaged  density  of  state 
with  energy  level  can  be  due  to  the  existence  of  indium-composition-fluctuation  nanostructures, 
which  is  caused  by  the  spinodal  decomposition  process,  in  the  sample.  ©  2005  American  Institute 
of  Physics.  [DOI:  10.1 063/1 . 1 847705] 


I.  INTRODUCTION 

Because  of  the  large  lattice  mismatch  between  GaN  and 
InN,  two  important  phenomena  in  InGaN  compounds  have 
been  widely  observed,  including  indium-composition  fluc¬ 
tuation  (or  clustering)1"4  and  strain-induced  piezoelectric 
field  near  an  interface  of  a  heterostructure.5"9  The  former 
leads  to  carrier  (exciton)  localization  and  hence  radiative  ef¬ 
ficiency  improvement.  The  latter  results  in  quantum-confined 
Stark  effect  (QCSE)  and  hence  carrier  separation  in  a  carrier- 
confined  structure.  Both  phenomena  are  important  in  inter¬ 
preting  the  photon  emission  mechanisms  in  such  a  com¬ 
pound.  In  particular,  these  two  phenomena  interplay  in  an 
InGaN/GaN  quantum  well  (QW)  structure,  leading  to  com¬ 
plicated  optical  behaviors.  Normally,  when  the  hetero¬ 
structure-induced  strain  is  relaxed  through  a  certain  process, 
spinodal  decomposition  can  be  easily  induced  such  that  in¬ 
dium  composition  fluctuations  and  indium-rich  clusters  can 
be  formed  in  such  a  sample.10,11  On  the  other  hand,  if  the 
strain  is  preserved,  the  solid  miscibility  between  InN  and 
GaN  can  be  increased  such  that  relatively  more  uniform  In¬ 
GaN  temaiy  compounds  can  exist.  In  this  situation,  higher 
QW  quality  can  be  obtained. 

With  the  effects  of  carrier  localization  and  QCSE  mixed 
in  a  QW  sample,  it  is  usually  difficult  to  identify  either  effect 
clearly.  However,  in  an  InGaN  thin-film  sample  of  large 
enough  thickness  (larger  than  the  critical  thickness),  the 
strain  effect  can  be  relaxed  in  most  part  of  the  sample.  In  this 
situation,  we  can  exclude  the  QCSE  and  observe  the  effect  in 
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the  optical  property  of  potential  fluctuation  solely.  To  under¬ 
stand  the  photon  emission  mechanisms  of  such  a  material, 
ultrafast  carrier  dynamics  observed  with  femtosecond  pump- 
probe  experiments  have  been  proven  an  effective  tool.12,13 
Carrier  flow  information  can  help  us  in  understanding  the 
nanostructures,  which  provide  us  with  the  clues  for  light- 
emission  origins.  In  this  paper,  we  report  the  results  of  ul¬ 
trafast  carrier  dynamics  in  an  InGaN  thin  film  with  the  thick¬ 
ness  larger  than  the  critical  thickness,  particularly  near  the 
energy  levels  of  free-carrier  states.  Temperature-,  pump- 
photon-energy-,  and  pump-intensity-dependent  femtosecond 
degenerate  pump-probe  differential  transmission  traces  show 
significant  variations  of  important  implications.  Such  varia¬ 
tions  are  attributed  to  the  variation  of  the  space-averaged 
density  of  state,  which  is  due  to  the  indium-composition 
fluctuation  in  the  sample.  This  paper  is  organized  as  follows: 
Descriptions  of  the  sample  structure  and  experimental  meth¬ 
ods  are  presented  in  Sec.  II.  The  pump-probe  experimental 
results  are  discussed  in  Sec.  III.  More  discussions  are  given 
in  Sec.  IV.  Finally,  conclusions  are  drawn  in  Sec.  V. 

II.  SAMPLE  PREPARATION  AND  EXPERIMENTAL 
METHODS 

The  InGaN  thin-film  sample  was  grown  on  c-plane  sap¬ 
phire  with  metal-organic  chemical-vapor  deposition 
(MOCVD).  After  the  800-nm  GaN  buffer  layer,  an  800-nm 
InGaN  thin  film  with  silicon  doping  of  5  X  1018  cm"3  in  con¬ 
centration  was  grown  at  800  °C.  The  average  indium  content 
was  estimated  to  be  20%.  Because  the  critical  thickness  of 
InGaN  under  the  current  conditions  is  only  about  50  nm,14 
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FIG.  1.  A  HRTEM  image  of  the  sample. 


the  results  of  transmission  measurement  mainly  originate 
from  the  strain-relaxed  portion  of  InGaN.  However,  compo¬ 
sition  fluctuations  of  roughly  10  nm  in  scale  size  were  ob¬ 
served  in  the  sample  with  high-resolution  transmission  elec¬ 
tron  microscopy  (HRTEM),  as  shown  in  Fig.  1.  The  bright 
regions  in  the  HRTEM  image  represent  the  indium-rich  dis¬ 
tributions.  The  HRTEM  investigation  was  performed  using  a 
Philips  TECNAI  F20  field-emission  electron  microscope 
with  an  accelerating  voltage  of  200  kV.  The  high-resolution 
micrograph  was  taken  at  the  Scherzer  defocus.  The  sample 
was  viewed  along  a11”20  zone  axis  for  a  bright-field  image. 

In  the  optical  experiment,  the  degenerate  pump  and 
probe  sources  originated  from  the  second-harmonic  of  a 
mode-locked  Ti: sapphire  laser  with  a  beta  barium  borate 
(BBO)  crystal.  The  second-harmonic  pulse  width  was  about 
150  fs.  The  tuning  range  was  from  2.7  through  3.4  eV.  The 
average  power  for  the  pump  ranged  from  4  through  42  mW. 
That  for  the  probe  was  fixed  at  0.5  mW.  In  the  temperature- 
dependent  measurements,  the  pump  average  power  was  fixed 
at  9  mW.  The  pump-probe  light  spot  on  the  sample  surface 
was  about  0.1  mm  in  diameter. 

IIS.  PUMP-PROBE  EXPERIMENTAL  RESULTS 

Figure  2  shows  the  temperature-dependent  variation  of 
photoluminescence  (PL)  spectra  and  the  spectral  peak  posi¬ 
tions  as  functions  of  temperature.  Two  PL  spectral  peaks  are 
identified  below  200  K.  An  S-shaped  variation  for  the  low- 
energy  peak  can  be  clearly  seen.  Such  a  variation  has  been 
regarded  as  one  of  the  optical  features  of  indium-rich  clusters 
in  the  InGaN  thin  film.3,15  The  low-energy  peak  corresponds 
to  the  localized  states.  The  high-energy  peak  is  attributed 
to  the  free-carrier  states,  corresponding  to  the  background 
InGaN  compound,  on  which  clusters  are  distributed.  How¬ 
ever,  it  is  noted  that  the  background 
InGaN  compound  also  consists  of  potential  fluctuations  of 
shallower  distributions.  The  merge  of  the  two  PL  peaks 
above  200  K  is  attributed  to  carrier  liquidation  among  the 
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FIG.  2.  Temperature-dependent  PL  spectra  and  the  calibrated  spectral  peak 
positions  for  the  high  and  low  energy  peaks. 

localized  states  and  free-carrier  states.  The  decreasing  trend 
of  the  high-energy  peak  is  mainly  due  to  the  band-gap 
shrinkage  of  phonon  effect.  The  concerned  photon  energy 
range  in  our  pump-probe  experiments  is  close  to  the  high- 
energy  peak  of  PL  spectrum.  In  other  words,  we  are  inter¬ 
ested  in  the  carrier  dynamics  at  the  energy  levels  around 
free-carrier  states. 

Figure  3  shows  the  differential  transmission  (A  777) 
traces  of  various  temperatures  with  the  pump  photon  energy 
at  3.2  eV.  The  time-resolved  traces  show  quite  different  be¬ 
haviors  as  temperature  increases  from  10  through  300  K. 
Below  30  K,  after  the  strong  pump-induced  peak,  the  probe 
intensity  decays  with  a  double-exponential  profile  within  15 
ps.  The  pump-induced  peak  at  30  K  is  particularly  strong  (up 
to  14%),  indicating  the  strong  band-filling  effect  at  this  tem¬ 
perature.  However,  between  60  and  180  K,  the  pump- 
induced  probe  intensity  shows  a  sharp  decrease  first  and  then 
an  increase  to  a  certain  level.  This  level  is  lower  than  that  of 
negative  delay  (without  pump)  in  either  trace  of  90  and  1 20 
K.  In  this  temperature  range,  the  A  77  T  variation  is  the  weak¬ 
est  among  the  traces  in  Fig.  3.  In  the  temperature  range  be¬ 
tween  60  and  1 80  K,  the  probe  intensity  decreases  and  then 
increases  from  the  peak  of  the  aforementioned  level.  In  some 
of  these  traces,  the  levels  around  15  ps  are  higher  than  those 
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FIG.  3.  Differential  transmission  traces  of  different  sample  temperatures 
with  pump  photon  energy  at  3.2  eV.  The  pump  pulse  energy  density  is 
3.26  juJ/cm2.  Two  zero  delay  times  for  different  traces  are  shown  with  the 
vertical  dashed  lines. 
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FIG.  4.  Differential  transmission  traces  of  different  sample  temperatures 
with  pump  photon  energy  at  3.14  eV.  The  pump  pulse  energy  density  is 
3.26  joJ/cm2. 

of  negative  delay.  As  temperature  increases  beyond  180  K, 
the  behavior  of  double-exponential  decay  from  the  pump- 
induced  peak  reappears.  The  peak  level  reaches  another 
maximum  (about  9%)  around  210  K.  In  the  traces  of  30  and 
210  K,  quasiperiodical  oscillations  can  be  observed.  Such 
oscillations  correspond  to  a  frequency  of  0.54  THz.  They  are 
attributed  to  the  interaction  of  coherent  acoustic  phonon.16"18 
Such  coherent  phonon  interaction  is  expected  at  all  tempera¬ 
tures;  however,  we  can  observe  it  only  at  those  of  significant 
band-filling  effects.  The  acoustic  phonon  is  created  at  the 
incidence  of  the  pump  laser.  In  a  multiple  QW  sample,  it 
travels  through  the  sample  and  interacts  with  the  carriers  in 
the  wells.  Therefore,  the  oscillation  period  of  the  pump- 
probe  trace  corresponds  to  the  period  of  the  QW.18  In  our 
case,  by  assuming  that  the  acoustic  speed  in  the  InGaN  thin 
film  is  6000  m/s,  the  observed  oscillations  of  around  0.54 
THz  correspond  to  a  periodical  structure  of  1 1  nm  in  period. 
In  the  used  sample,  although  composition  fluctuations  are 
observed  in  the  HRTEM  image  (see  Fig.  1),  we  have  no 
strong  evidence  for  this  periodical  structure.  However,  there 
is  no  reason  to  exclude  such  a  possibility  either.  This  scale 
size  (around  1 1  nm)  may  represent  a  characteristic  dimension 
of  the  indium-rich  clusters  or  composition  fluctuations.  Nev¬ 
ertheless,  it  deserves  further  investigation. 

Figure  4  shows  the  AT/T  traces  at  various  temperatures 
with  the  pump  photon  energy  at  3.14  eV.  Here,  one  can  ob¬ 
serve  the  behavior  of  double-exponential  decay  after  the 
pump-induced  peak  in  the  temperature  range  below  180  K. 
Above  this  temperature,  the  behavior  of  the  probe  intensity 
decrease  at  the  pump  instant  can  be  seen.  The  trace  of  210  K 
is  similar  to  that  of  120  K  in  Fig.  3.  However,  above  this 
temperature,  the  small  feature  rise  after  the  dip  disappears. 
Here,  a  double-exponential  rising  behavior  can  be  observed. 

To  summarize  the  data  in  Figs.  3  and  4,  and  the  similar 
results  of  other  pump  photon  energies,  the  photoinduced 
change  of  optical  density  (AOD)  versus  temperature  at  three 
pump-probe  photon  energies  is  shown  in  Fig.  5. 19  The  optical 
density  is  defined  as  the  ratio  of  unpumped  transmission  over 
the  pumped  transmission  in  the  log  scale  at  zero  pump-probe 
delay,  i.e.,  the  temporal  positions  indicated  with  the  vertical 
dashed  lines  in  Figs.  3  and  4.  In  Fig.  5,  variations  of  AOD  of 


FIG.  5.  AOD  vs  temperature  of  three  pump  photon  energies  and  AOD  vs 
photon  energy  at  room  temperature. 

different  amplitudes  for  the  three  photon  energies  can  be 
observed.  In  the  inset  of  Fig.  5,  AOD  data  versus  photon 
energy  at  room  temperature  are  plotted.  Here,  a  quasiperiodi¬ 
cal  variation  with  a  period  around  120  meV  can  also  be  seen. 

Figure  6  shows  the  A  TiT  traces  of  various  pump  pulse 
energy  densities  (equivalent  to  pump  intensities)  at  room 
temperature  with  the  pump  photon  energy  at  3.1  eV.  One  can 
see  that  as  the  pump  pulse  energy  density  increases  from 
2.89  fill  cm2  (corresponding  to  1 .94  X  1017  cm-3  photoge¬ 
nerated  carrier  density)  to  16.53  /LtJ/cm2,  the  AT/ T  profile 
changes  from  the  pattern  of  single  hump  into  that  of  a  dip. 

IV.  DISCUSSIONS 

A  77  T  represents  the  transient  absorption  alteration  of  the 
probe  upon  the  pump  excitation.  A  trace  with  a  rise,  followed 
by  a  double-exponential  decay  (including  the  faster  local  car¬ 
rier  thermalization  plus  carrier  relaxation  into  lower  energy 
levels  and  the  slower  carrier  diffusion  or  global  carrier 
thermalization17),  is  normally  interpreted  as  the  band-filling 
effect,  which  diminishes  with  the  carrier  flow  and  recombi¬ 
nation.  In  such  a  process,  band-gap  renormalization  and  pho¬ 
non  effect  can  also  occur.  With  a  high  localized-carrier  den¬ 
sity  upon  the  pump  excitation,  band-gap  renormalization 
leads  to  a  reduced  band  gap.  This  process  combines  with  the 
phonon  effect  for  the  resultant  effective  band-gap  alteration. 


FIG.  6.  Differential  transmission  traces  of  different  pump  photon  energy 
densities  at  room  temperature  with  the  pump  photon  energy  at  3.1  eV. 
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When  the  band  gap  varies,  the  pump-probe  photon  corre¬ 
sponds  to  the  energy  levels  of  different  space-averaged  den¬ 
sities  of  state,  which  determine  the  band-filling  processes  and 
hence  the  A  TIT  behaviors.  Therefore,  the  variations  of  the 
A  TiT  behavior,  in  changing  either  temperature,  pump  photon 
energy,  or  pump  intensity,  should  follow  the  variation  pattern 
of  the  space-averaged  density  of  state  along  the  energy  level. 
Here,  the  space-averaged  density  of  state  means  the  average 
energy-level  density  over  a  spatial  region  of  composition 
fluctuations.  In  this  region,  because  of  the  spatial  potential 
variation,  random  or  organized,  space-averaged  density  of 
state  may  show  strong  dependence  on  the  energy  level.  A 
quasiperiodical  variation  can  exist  in  such  a  structure  formed 
through  spinodal  decomposition.  As  shown  in  Figs.  3  and  4, 
when  temperature  increases,  the  band  gap  shrinks  and  the 
pump-probe  photon  energy  moves  to  the  higher  energy  level 
with  respect  to  the  band  edge.  In  this  situation,  the  pump 
may  experience  different  densities  of  state  and  hence  differ¬ 
ent  absorption  coefficients  such  that  the  A  77  T  shows  differ¬ 
ent  magnitudes  with  increasing  temperature.  Also,  when  we 
change  the  pump-probe  photon  energy,  as  shown  in  Fig.  5,  it 
is  expected  to  observe  different  optical  densities  at  a  fixed 
temperature.  Meanwhile,  with  localized  band-gap  renormal¬ 
ization,  the  induced  band-gap  shrinkage  and  potential- 
fluctuation  reform  may  change  the  corresponding  space- 
averaged  density  of  state  and  hence  the  pump-probe  behavior 
when  the  pump  intensity  is  increased,  as  shown  in  Fig.  6. 

When  the  corresponding  density  of  state  is  high,  a  sig¬ 
nificant  absorption  leads  to  a  strong  band-filling  effect,  re¬ 
sulting  in  the  high  peaks  around  the  zero  pump-probe  delay 
and  the  two-component  decays.  On  the  other  hand,  when  the 
corresponding  density  of  state  is  low,  the  band-filling  effect 
is  weak  and  hence  two-photon  absorption  (TPA)  and  free- 
carrier  absorption  (FCA)  become  observable.20  The  sharp  de¬ 
creases  or  dips  in  some  of  the  A 77  T  traces  in  Figs.  3,  4,  and 
6  are  attributed  to  TPA.  In  these  situations,  the  decays  of  the 
peaks  right  after  the  dips  or  the  elongations  of  the  dips  are,  at 
least  partly,  due  to  FCA.  It  is  noted  that  TPA  and  FCA  are 
expected  to  occur  in  any  situation,  either  high  or  low  density 
of  state.  However,  they  are  prominent  only  when  the  band¬ 
filling  effect  is  weak. 

V.  CONCLUSIONS 

In  summary,  we  have  performed  femtosecond  degenerate 
pump-probe  experiments  on  an  InGaN  thin  film.  The  ob¬ 
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served  temperature-,  pump-photon-energy-,  and  pump- 
intensity-dependent  variations  of  ultrafast  carrier  dynamics 
manifested  the  variation  of  space-averaged  density  of  state 
along  the  energy  level  in  the  sample,  in  which  potential  fluc¬ 
tuations  existed  due  to  indium-composition  fluctuations.  The 
carrier  dynamics  was  controlled  by  the  shift  of  effective  band 
gap  and  hence  the  behavior  of  band  filling.  The  shift  of  ef¬ 
fective  band  gap  was  determined  by  the  combined  effect  of 
band-gap  renormalization  and  phonon  effect.  TPA  and  FCA 
became  observable  when  the  density  of  state  was  low  and 
hence  the  band-filling  effect  was  weak.  For  the  future  work, 
the  detailed  variation  of  the  space-averaged  density  of  state 
in  an  InGaN  thin  film  deserves  further  investigation. 
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I  1.  Introduction 

It  has  been  shown  that  different  doping  condi¬ 
tions,  including  different  doping  layers,  in  InGaN/ 
5  GaN  quantum  well  (QW)  structures  led  to 
different  optical  characteristics.  Such  effects  of 
7  silicon  doping  have  been  widely  studied.  Several 
models  were  proposed  for  interpreting  the  optical 
9  and  material  characteristics  of  InGaN/GaN  QWs 
with  silicon  doping.  The  most  commonly  used 
1 1  model  is  the  carrier  screening  effect  of  the  strain- 
induced  piezoelectric  field  and  hence  the  reduction 
13  of  the  quantum-confined  Stark  effect  (QCSE). 

This  effect  leads  to  the  reduced  Stokes  shift  [1-3], 
15  the  decrease  of  radiative  recombination  lifetime 
[2-4,10],  the  blue  shift  of  PL  spectrum  [3,5],  and 
17  the  enhancement  of  photon  emission  efficiency 
[1,3,6].  Also,  material  analyses  have  led  to  the 
19  conclusions  of  growth  mode  change  [6],  nanos¬ 
tructure  alteration  [7],  formation  of  abrupt  QW 
21  interfaces  [2,8],  strain  relaxation  [7,9],  and  higher 
potential  uniformity  in  QWs  [3].  Meanwhile,  it  has 
23  been  shown  that  different  barrier-doping  concen¬ 
trations  in  InGaN/GaN  QW  structures  resulted  in 
25  different  optical  characteristics  [11].  These  results 
are  not  necessarily  mutually  consistent,  particu- 
27  larly  in  the  material  nanostructures  of  such  silicon- 
doped  QW  samples. 

29  In  particular,  due  to  the  large  lattice  mismatch 
between  InN  and  GaN,  indium-rich  clusters  are 
31  usually  formed  in  such  a  QW  structure, ftjirough 
spinodal  decomposition  and  other  effects  [12-14]. 
33  In  some  situations,  quantum-dot-like  nanostruc¬ 
tures  can  be  observed  [14 — 16].  With  different 
35  silicon  doping  conditions,  the  clustering  structures 
are  expected  different.  Recently,  the  techniques  of 
37  high-resolution  transmission  electron  microscopy 
(HRTEM)  have  been  widely  used  for  understand- 
39  ing  the  properties  of  threading  dislocations  [17,18] 
and  nano-structures  in  GaN  and  InGaN/GaN 
41  QW  samples.  In  particular,  the  existence  of 
indium-rich  clusters  in  InGaN/GaN  QW  struc- 
43  tures  has  been  an  interesting  issue  of  discussion 
[19-22].  Although  some  researchers  claimed  that 
45  the  observations  of  clusters  were  due  to  the 

artifacts  [19,20],  created  through  electron  beam 
47  exposure  during  HRTEM  measurements,  the 

existence  of  clusters  really  relies  on  the  conditions 


of  sample  preparation.  Under  different  conditions  49 
of  crystal  growth,  clusters  may  or  may  not  exist. 

Also,  the  cluster  characteristics  depend  on  those  51 
conditions. 

In  this  paper,  we  use  the  strain-state  analysis  53 
(SSA)  method  [23]  to  study  the  nanostructures  of 
InGaN/GaN  QWs  of  different  doping  conditions  55 
(un-doped,  well-doped,  and  barrier-doped)  and 
different  average  indium  contents.  It  is  found  that  57 
the  clustering  phenomenon  is  the  strongest  in  the 
barrier-doped  samples  and  the  weakest  in  the  un-  59 
doped  samples.  Among  the  samples  of  different 
average  indium  contents,  clustering  behavior  is  61 
relatively  stronger  in  the  high-indium  samples. 

Based  on  the  SSA  results,  the  detailed  nanostruc-  63 
tures  are  analyzed  and  the  relative  indium  contents 
can  be  confirmed.  Besides,  optical  properties  based  65 
on  photoluminescence  (PL)  measurements  are 
provided  to  show  the  consistent  trends  in  the  67 
comparison  of  different  doping  conditions.  This 
paper  is  organized  as  follows:  in  Section  2,  the  69 
sample  structures,  growth  conditions,  and  mea¬ 
surement  procedures  are  reported.  The  SSA  results  71 
of  various  samples  are  discussed  in  Section  3. 

Then,  the  optical  characterization  results  are  73 
presented  in  Section  4.  Certain  discussions  are 
given  in  Section  5.  Finally,  conclusions  are  drawn  75 
in  Section  6. 

77 

2.  Sample  descriptions  and  experimental  procedures  79 

The  six  InGaN/GaN  QW  samples  of  different  81 
doping  conditions  and  different  average  indium 
contents  were  prepared  with  MOCVD  growth  on  83 
sapphire  (0001)  substrate.  They  all  consisted  of 
five  QW  periods,  with  2.5 nm  in  well  width  and  85 
7.5  nm  in  barrier  width,  on  top  of  an  un-doped 
GaN  buffer  layer  of  1.52  pm  in  thickness.  The  87 
growth  temperatures  were  1 100  and  800  °C  for  the 
GaN  barriers  and  InGaN  wells,  respectively.  The  89 
doping  concentration  of  silicon  was  5  x  1018cm-3, 
either  in  barriers  or  wells.  They  are  classified  into  91 
two  groups  of  low  and  high  average  indium 
contents.  The  low-indium  set  of  samples  has  an  93 
estimated  content  of  15%.  The  high-indium  set  has 
indium  content  of  25%  in  estimation.  In  the  low-  95 
indium  set,  the  samples  of  un-doped,  well-doped 
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1  and  barrier-doped  conditions  were  denoted  with 
samples  LU,  LW,  and  LB,  respectively.  The 
3  counterparts  of  the  high-indium  set  are  assigned 
as  samples  HU,  HW,  and  HB,  respectively. 

5  The  HRTEM  investigations  (for  SSA)  were 
performed  using  a  Philips  Tecnai  F30  field- 
7  emission  electron  microscope  using  an  accelerating 
voltage  of  300  kV  and  a  probe-forming  lens  of 
9  Cs  =  1.2  mm.  All  the  high-resolution  micrographs 
were  taken  with  the  two-electron-beam  interfer- 
1 1  ence.  The  current  density  focused  onto  the  samples 
was  estimated  to  be  lower  than  16  A/cm2.  With  the 
13  two-beam  interference  data,  we  could  conduct  the 
SSA  for  composition  distribution  images  [23].  In 
15  the  HRTEM  operation,  to  avoid  the  possible 
influence  of  electron  illumination  on  the  recorded 
17  results  [19,20],  the  two-beam  interference  images 
were  always  taken  within  5-10  s  after  the  electron 
19  beam  was  focused  onto  the  target  spots.  Our 
experiments  showed  that  the  processed  SSA 


21  images  were  not  significantly  changed  within  30s  — Jl_ 69 

of  the  illumination  of  focused  electron  beams.  To  l*°°  ,*01  1,02  ,#03  1,04  I#05  107  >l’°* 

23  obtain  the  SSA  images,  the  [0001]  lattice  size  Xfc^%)  0  6  13  I$  24  31  36  42  >4$  7] 

distributions  were  first  calibrated  from  the  lattice  Fig,  j.  Typical  hrtem  lattice  fringe  image  (a)  and  the  ssa 

25  fringe  images.  Then,  based  on  Vegard’s  law  for  the  ^:4  imag^)  after  software  process  of  samples  lu.  73 

relation  between  the  lattice  size  and  indium 

27  composition,  composition  images  could  be  ob-  75 


tained.  These  compositions  were  estimated  based  ,  parameter  ratio  dj/do snn  (^GaN  and  dj  for  the 
29  on  the  assumption  of  a  specimen  thickness:  larger  lattice  sizes  of  GaN  and  InGaN,  respectively)  and  77 

than  50  nm.  PL  measurements  were  carried out  the  indium  composition  xin  (%)  is  shown  in  the 

3 1  with  the  325  nm  line  of  a  35  mW  He-Gd  iaser.  The  legend  at  the  bottom  of  the  figure.  It  is  noted  that  79 

samples  were  placed  in  a  cryostat  for  temperature-  their  relationship  is  not  linear,  particularly  when 

33  dependent  measurements  ranging  from  10  to  the  indium  composition  is  high.  In  each  SSA  81 

300  K.  ,.,:K  image,  line  scan  profiles  of  indium  composition 

35  along  the  indicated  white  lines  are  plotted.  The  83 

scales  for  the  profiles  in  the  two  dimensions  are  the 

37  3.  Images  of  SSA  same.  In  each  SSA  image,  a  QW  surrounded  by  85 

'  ■ ^  v  the  two  barriers  are  illustrated. 

39  Figs.  1-6  show  the  typical  SSA  results  of  Fig.  1(b)  shows  the  SSA  image  of  sample  LU.  87 
samples  LU,  LW,  LB,  HU,  HW,  and  HB,  One  can  see  that  the  QW  layer  (green  region)  is 

41  respectively.  In  each  figure,  (a)  shows  the  lattice  quite  uniform.  The  indium  clusters  (reddish  or  89 

fringe  image  of  two-beam  interference  along  the  dark  regions)  within  the  QW  layer  are  quite  small 

43  [000  1]  direction  and  (b)  shows  the  processed  SSA  and  shallow.  The  GaN  barriers  (deep  blue  or  91 

image.  One  can  see  that  the  lattice  fringe  image  bright  regions)  are  clear.  It  means  that  indium  is 

45  cannot  clearly  describe  the  nanostructures  of  well  confined  within  the  well.  From  the  horizontal  93 

materials.  The  colors  (gray  levels)  in  the  SSA  line-scan  profile,  one  can  observe  quite  a  weak 

47  image  represent  the  estimated  indium  composition  fluctuation  in  indium  composition  along  the  well  95 

distributions.  The  relationship  between  the  lattice  layer.  The  fluctuation  contrast  (the  difference 
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Fig.  2.  Typical  HRTEM  lattice  fringe  image  (a)  and  the  SSA 
image  (b)  after  software  process  of  samples  LW. 

between  the  maximum  and  the  minimum  Jn  the 
scan  range)  is  in  the  range  of  0.03-0.04  (18-24%) 
along  the  QW  layer.  In  Fig.  2(b)  for  sample  LW, 
we  can  observe  that  more  indium  clusters  are 
formed  in  the  QW  layer.  Also,  more  indium  atoms 
diffuse  from  the  well  into  barriers*  Here,  the 
fluctuations  of  the  line-scan  profiles  are  relatively 
stronger  in  comparing  with  those  of  sample  LU. 
The  fluctuation  contrast  is  now  in  the  range  of 
0.04-0.05  (22-28%).,  In  Fig.  3(b)  for  sajnple  LB, 
we  can  see  two  clear  indium- rich  clusters  in  the 
QW  (two  largef^edi^pots).  From  the  line-scan 
profiles,  one%n%e  that  the  indium  composition 
fluctuatioir  in  sample  LB  is  much  stronger  than 
those  of  thlSfeffier  two  samples.  The  fluctuation 
contrast  is  now  around  0.08  (44%). 

As  shown  in  Fig.  4(b)  for  sample  HU,  although 
the  QW  interface  is  blurred,  indium  is  basically 
confined  within  the  well.  Here,  within  the  well  a 
red  spot  of  indium  aggregation  can  be  observed. 
Similar  to  sample  LU,  one  can  also  observe  a 


d/d^  1.00  1.01  1.02  1.03  1.04  1.05  1 M  147  >U0S 
0  6  13  1$  34  31  36  42  >4$ 


HRTEM  lattice  fringe  image  (a)  and  the  SSA 
^imagiltQ  after  software  process  of  samples  LB. 

% 

relatively  weaker  fluctuation  in  indium  composi¬ 
tion  along  the  well  layer  in  the  line-scan  profile. 
The  fluctuation  contrast  is  about  0.04  (22%)  along 
the  QW  layer.  In  Fig.  5(b)  for  the  SSA  image  of 
sample  HW,  the  QW  is  not  as  well  shaped  as 
sample  HU.  The  indium  composition  fluctuations 
in  both  directions  are  relatively  stronger  in 
comparing  with  sample  HU.  In  particular,  more 
indium-aggregated  clusters  can  be  observed  within 
the  well  layer.  The  fluctuation  contrast  is  also 
around  0.04  (22%)  along  the  QW  layer.  As  shown 
in  Fig.  6(b),  the  SSA  image  of  sample  HB  shows 
quite  a  different  nanostructure  from  the  other  five 
samples.  Here,  the  QW  layer  becomes  unclear. 
Instead,  a  distribution  of  clusters  of  different  sizes 
and  shapes  exists.  The  indium  composition  fluc¬ 
tuation  is  much  stronger  than  those  of  samples 
HU  and  HW.  The  contrast  is  now  around  0.08 
(44%)  along  the  QW  layer. 

A  question  has  been  raised  about  the  possible 
influence  of  electron  beam  exposure  on  the 
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Fig.  4.  Typical  HRTEM  lattice  fringe  image  (a)  and  the  SSA 
image  (b)  after  software  process  of  samples  HU. 


Typical  HRTEM  lattice  fringe  image  (a)  and  the  SSA 
after  software  process  of  samples  HW. 


im¬ 


material  nanostructures  during  the  measurement 
of  transmission  electron  microscopy  [19,20]ypb 
have  tested  such  a  possibility  by  varying  "the 
electron  exposure  time.  The  SSA  images  of  sample 
LB  shown  in  Fig.  7(a)  and  (b)  ^g^^^Ted  under 
electron  beam  exposure  Jp^^iM|ian  10  s  and 
about  2  min,  respectivel>^FiM7i^)  is  different 
from  Fig.  3(b)  becai^Sr  m^mere  obtained  at 
different  locations  of||amp||  LB.  Between  Fig.  7(a) 
and  (b),  one  can  ^e  ®a^Slthough  the  details  are 
slightly  changed^the  basic  structures  are  the  same. 
In  other  words,  with  up  to  2-min  exposure,  we  can 
still  obtain  the  required  information  about  the 
nanostructures  of  the  samples.  Normally,  when  we 
operated  the  HRTEM  system,  it  took  less  than  10  s 
after  the  electron  beam  was  focused  onto  a  target 
spot  to  obtain  an  image.  In  our  research,  we  have 
done  our  best  to  reduce  the  electron  beam  effect. 
Hence,  the  SSA  images  in  this  research  can  truly 
represent  the  sample  nanostructures. 


In  Figs.  8  and  9,  we  show  the  average  line-scan 
profiles  in  the  [0001]  direction,  averaged  over  the 
image  range  along  the  [11-20]  direction,  plotted  as 
functions  of  distance  for  the  low-  and  high-indium 
samples,  respectively.  Two  to  three  SSA  images 
like  Figs.  1-6  at  different  locations  in  each  sample 
were  used  for  taking  the  average.  One  can  see  that 
generally  the  well  widths  of  the  high-indium 
samples  are  larger  than  those  of  the  low-indium 
ones.  One  can  also  observe  that  the  interfaces 
between  the  well  and  barriers  of  the  low-indium 
samples  are  generally  sharper  than  those  of  the 
high-indium  ones.  Meanwhile,  generally  the  in¬ 
dium  concentrations  at  the  peaks  of  the  high- 
indium  samples  are  higher  than  those  of  the  low- 
indium  ones.  In  addition,  in  the  low-indium 
samples,  more  indium  atoms  diffuse  from  the  well 
into  the  lower  barrier.  Indium  out-diffusion  in  the 
high-indium  samples  is  more  symmetric. 

The  average  indium  content  of  each  sample  can 
be  calibrated  by  integrating  the  average  line-scan 
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Fig.  8.  Average  indium  composition 


along  the  crystal 


growth  direction  for  samples  LU,  LW,  and  LB. 


Fig.  9.  Average  indium  composition  profiles  along  the  crystal 
growth  direction  for  samples  HU,  HW,  and  HB. 


town 

^  tiiaq 

fiTOOi 


results  in  Figs.  8  or  9  over  the  whole  range  and 
then  dividing  the  integrated  result  by  the  desig¬ 
nated  QW  width  (2.5  nm).  The  whole-range 
integration  implies  the  coverage  of  the  indium 
atoms,  which  diffuse  from  the  well  into  barriers. 
Such  results  of  average  indium  content  become 
more  accurate  if  more  SSA  images  like  Figs.  1-6 
are  used  for  average.  Although  the  current  results 
cannot  be  very  accurate,  they  can  still  be  used  for  a 
general  comparison.  The  calibrated  average  in¬ 
dium  contents  of  samples  LU,  LW,  and  LB  are 
14.3%,  17.1%,  and  19.2%,  respectively.  Those  of 
samples  HU,  HW,  and  HB  are  25.6%,  26.7%,  and 
26.9%,  respectively.  These  results  are  generally 
consistent  with  the  designated  indium  contents  of 
15%  and  25%  for  the  low-  and  high-indium 


t  .00  1.01  1.02  1.03  1.04  1.05  1.06  1.07  >1.0* 

0  6  13  1*  24  31  36  42  >48 

Fig.  6.  Typical  HRTEM  lattice  fringe  image  (a)  and  the  SSA 
image  (b)  after  software  process  of  samples  HB. 


[0001] 


[llOO]  [1150] 


d/d^  1.00  1.01  1.02  1.03  1.04  1  05  1.06  1.07  >1.00 
*■»  0%  6%  13%  1S%  24*  31*  36*  42*  >40* 

Fig.  7.  SSA  images  of  sample  LB  under  different  electron-beam 
exposure  times:  (a)  <10s  and  (b)  120  s. 
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samples,  respectively.  The  generally  higher  indium 
contents  from  calibration  for  the  barrier-  and  well- 
doped  samples  can  be  attributed  to  the  stronger 
indium  diffusion  into  barriers  in  these  two  kinds  of 
sample.  With  more  indium  atoms  in  the  barriers, 
the  calibration  for  indium  composition  in  proces¬ 
sing  the  SSA  images  can  be  overestimated. 
However,  the  possibility  of  really  higher  average 
indium  content  in  the  doped  samples  cannot  be 
ruled  out.  With  different  strain  conditions  (due  to 
different  doping  conditions),  the  incorporation 
efficiency  of  indium  atom  can  be  different  even 
the  flow  rate  is  the  same  during  growth. 


4.  Optical  characterization  results 

To  understand  the  corresponding  optical  prop¬ 
erties  of  these  samples,  basic  optical  characteriza¬ 
tions  were  performed.  Fig.  10  shows  the  PL 
spectral  peak  positions  as  functions  of  temperature 


Fig.  10.  PL  spectral  peak  variations  with  temperature  of  the  six 
samples. 


-i  — ,^,r.  i _ /  1  .  t  _ L 
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Fig.  11.  Integrated  PL  intensity  variations  with  temperature  of 
samples  LU,  LW,  and  LB. 

for  all  the  six  samples  with  the  filled  and  empty 
symbols  for  and  high-indium  samples, 

respectively.  Clear  S-shaped  variations  can  be  seen 
in  the  curves  for  samples  LU  and  HU.  On  the 
contrary;  the  S-shaped  behaviors  are  unclear  in  the 
well-  and  barrier-doped  samples.  The  PL  peak 
positions  of  the  un-doped  and  well-doped  samples 
eari^ither  group  of  sample  are  quite  close, 
ipart^iilarly  in  the  low-temperature  range.  The 
significant  blue  shifts  of  the  PL  peaks  of  samples 
LB  and  HB,  with  respect  to  those  of  the  un-doped 
I  and  well-doped  samples  show  a  clear  feature  of 
barrier  doping. 

Fig.  11  shows  the  integrated  PL  intensities  as 
functions  of  temperature  of  samples  LU,  LW,  and 
LB.  Fig.  12  shows  the  counterparts  of  samples 
HU,  HW,  and  HW.  The  thermal  quenching 
behavior  in  integrated  luminescence  intensity 
represents  the  radiative  efficiency.  From  Figs.  11 
and  12,  one  can  see  that  the  barrier-doped  samples 
always  have  the  highest  radiative  efficiency, 
followed  by  the  well-doped  samples  and  then  the 
un-doped  samples.  Hence,  in  our  samples,  silicon 
doping,  particularly  doping  in  barriers,  can  lead  to 
a  higher  efficiency  of  spontaneous  emission. 


5.  Discussions 

The  SSA  images  show  that  barrier  doping 
results  in  stronger  indium  clustering.  This  result 
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t  .  Fig.  12.  Integrated  PL  intensity  variations  with  temperature  of 
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may  be  caused  by  the  significant  relaxation  of 
19  hetero-structure-induced  strain  (the  strain  in  QW 
due  to  the  lattice  mismatch  between  wells  and 
21  barriers)  in  the  barrier-doped  samples,  particularly 
in  sample  HB,  through  certain  mechanisms.  In 
23  sample  HB,  with  the  growth  temperature  at  800  °C 
and  nominal  indium  content  about  25%,  the 
25  sample  condition  falls  into  the  regime  of  spinodal 
decomposition  (the  unstable  regime)  in  the  T-x 
27  phase  diagram  of  ternary  In^Gai-^N  compounds 
for  relaxed  and  strained  layers  with  the  interface 
29  orientation  perpendicular  to  the  hexagonal  axis  of 
the  crystal  [12,24],  In  this  regime,  spinodal 
31  decomposition  and  phase  separation  can  be 
initiated  by  a  small  local  compositional  fluctua- 
33  tion.  Hence,  strongly  clustering  behaviors  can  be 
observed  in  sample  HB.  The  mechanisms  for 
35  relaxing  the  hetero-structure-induced  strain  can 
be  the  generation  of  local  dislocations,  which  is 
37  caused  by  the  enhanced  lattice  mismatch  between 
a  well  and  the  surrounding  barriers  because  of  the 
39  substitutions  of  a  high-concentration  of  larger 
gallium  atoms  (atomic  weight  as  31)  by  smaller 
41  silicon  atoms  (atomic  weight  as  14)  in  the  barriers. 

The  HRTEM  images  in  the  similar  samples  of  our 
43  and  other’s  previous  research  have  shown  the 
existence  of  local  dislocations  in  an  InGaN/GaN 
45  QW  sample  of  a  high  nominal  indium  content 
[25,26],  The  other  mechanism  for  strain  relaxation 
47  in  the  barrier-doped  sample  is  the  existence  of 
plenty  carriers  in  the  barriers  and  wells  (carriers 


flow  into  wells)  due  to  the  heavy  silicon  doping  in  49 
the  barriers.  The  carriers  screen  the  piezoelectric 
fields  and  affect  the  strain  energy  distributions.  51 
However,  it  is  believed  that  the  major  mechanism 
for  strain  relaxation  is  the  formation  of  local  53 
dislocations.  The  discussions  described  above  fit 
well  the  conditions  of  sample  HB.  In  sample  LB,  55 
because  the  nominal  indium  content  is  lower 
(about  15%),  the  hetero-structure-induced  strain  57 
is  not  as  strong  as  that  in  sample  HB.  Therefore, 
the  condition  may  fall  into  a  region  between  the  59 
unstable  and  metastable  regimes  in  the  T-x  phase 
diagram.  In  this  situation,  spinodal  decomposition  61 
and  phase  separation  processes  can  still  be 
induced.  However,  these  processes  in  the  meta-  63 
stable  regime  take  a  longer  time  than  that  in  the 
unstable  regime,  leading  to  a  relatively  weaker  65 
clustering  behavior. 

In  the  well-doped  samples,  because  of  the  67 
substitutions  of  larger  gallium  or  indium  atoms 
by  smaller  silicon  atoms  (lattice  mismatch  re-  69 
duced)  and  the  effect  of  carrier  screening,  the 
hetero-structure-induced  strain  has  been  partly  71 

relaxed  such  that  the  conditions  may  belong  to  the 
meta-stable  regime  in  the  T-x  phase  diagram.  73 
Hence,  certain  clusters  were  generated  although 
the  clustering  behaviors  are  not  as  strong  as  those  75 
in  the  barrier-doped  samples.  In  the  un-doped 
samples,  the  strong  hetero-structure-induced  77 

strains  are  preserved  such  that  the  conditions  are 
essentially  located  in  the  stable  regime.  In  this  79 
situation,  the  solubility  of  InN  in  GaN  is  higher 
such  that  the  clustering  behaviors  are  the  weakest  81 
among  the  three  kinds  of  sample.  It  is  noted  that 
with  a  higher  nominal  indium  content,  the  condi-  83 
tion  in  the  phase  diagram  can  be  closer  to  the 
meta-stable  or  even  spinodal  decomposition  re-  85 
gime  and  hence  the  clustering  behavior  can  be 
stronger.  87 

From  the  SSA  results  above,  we  can  build  a 
model  for  potential  variation  along  a  QW  layer.  In  89 
the  situations  of  samples  HU  and  LU,  potential 
fluctuations  are  weak.  A  potential  minimum  is  91 
surrounded  by  secondary  minima  of  shallow 
barriers  in  between.  In  the  cases  of  samples  HW  93 
and  LW,  potential  variations  are  relatively  stron¬ 
ger.  In  this  situation,  a  potential  minimum  is  also  95 
surrounded  by  secondary  minima.  However,  the 
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1  barriers  between  the  minima  are  now  relatively 
higher  such  that  certain  energy  is  required  for 
3  carrier  transport  between  the  minima.  Then,  for 
samples  HB  and  LB,  sharp  minima  with  few 
5  secondary  dips  are  assumed.  Therefore,  carriers 
can  relax  down  to  the  absolute  potential  minima 
7  without  a  cascading  process  through  the  secondary 
minima.  The  difference  between  the  well-doped 
9  and  un-doped  samples  is  that  more  energy  is 
required  for  overcoming  the  barriers  in  samples 
1 1  LW  and  HW. 

The  relative  larger  PL  photon  energies  in  the 
13  barrier-doped  samples  can  be  due  to  the  reduction 
of  piezoelectric  fields  that  originates  from  strain 
15  relaxation  and  carrier  screening.  Also,  it  can  be 
due  to  the  increase  of  quantum  confinement, 
17  which  results  from  the  formation  of  nano-scale 
clusters.  In  the  well-doped  samples,  more  carriers 
19  in  the  well  layers  should  result  in  significant 
increases  of  PL  photon  energies  due  to  strong 
21  carrier  screening.  However,  the  PL  photon  en¬ 
ergies  of  the  well-doped  samples  are  quite  close  to 
23  those  of  the  un-doped  samples.  Such  a  result  can 
be  attributed  to  the  ineffective  carrier  screening 
25  effect.  It  can  also  be  interpreted  as  the  result  of  the 
counteraction  between  band  gap  renormalization 
27  (red-shift  trend)  and  piezoelectric  field  reduction 
(blue-shift  trend).  As  shown  in  Fig.  10,  the 
29  temperature-dependent  PL  spectral  peak,  energies 
reveal  clear  S-shaped  behaviors  in  samples  LU  and 
31  HU.  Such  an  S-shaped  behavior^,  is^^jjsually 
attributed  to  carrier  localization  [27,%8]  and/or 
33  QCSE  [29,30] .  Because  the  S-shaped  behavior  is 
unclear  in  either  well-doped  or  barrier-doped 
35  samples,  whose  strains  are  supposed  to  be  partially 
or  fully  relaxed,  the  major  mechanism  for  the  S- 
37  shaped  behavior  in  our  samples  must  be  the  strain 
distribution  or  QCSE.  Without  significant  relaxa- 
39  tion  of  strains  in  the  un-doped  samples,  the  S- 
shaped  behaviors  are  quite  prominent.  The  sig- 
41  nificant  enhancement  of  radiative  recombination 
efficiency  in  the  barrier-doped  samples  can  be 
43  attributed  to  the  strong  carrier  localization  (loca¬ 
lized  from  non-radiative  recombination  centers) 
45  and  the  reduced  QCSE  (improved  wave  function 
overlap  between  electron  and  hole). 

47 


6.  Conclusions  49 

In  summary,  we  have  compared  the  results  of  51 
SSA  and  PL  of  six  InGaN/GaN  QW  samples  with 
un-doped,  well-doped,  and  barrier-doped  struc-  53 
tures.  Based  on  the  SSA  images,  a  strain  relaxation 
model  was  built  for  describing  the  nanostructure  55 


differences  between  the  three  sets  of  sample  of 
different  doping  conditions.  In  the  barrier-doped  57 
samples,  the  hetero-structure-induced  strains  were 
fully  relaxed  such  that  spinodal  decomposition  59 
was  effectively  induced.  Therefore,  strongly  clus¬ 
tering  nanostructures  were  observed,  In  the  well-  61 
doped  samples,  strains  were  partially  relaxed  and 
the  spinodal  decomposition  process  could  be  63 
slightly  induced.  Hence,  weaker  composition 
fluctuations  were  observed.  Then,  in  the  un-doped  65 
samples,  the  un-relaxed  strains  resulted  in  higher 
miscibility  between  InN  and  GaN,  leading  to  the  67 
relatively  more  uniform  composition  distributions. 
Between  the  low-  and  high-indium  samples,  higher  69 
indium  content  led  to  a  stronger  clustering 
behavior.  The  strain  relaxations  in  the  well-doped  71 
and  barrier-doped  samples  resulted  in  their  unclear 
S-shaped  behaviors  of  PL  spectral  peaks.  The  73 
enhanced  carrier  localization  and  reduced  QCSE 
in  the  barrier-doped  samples  were  responsible  for  75 
their  significant  increases  of  radiative  efficiency. 
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